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ARTICLE INFO ABSTRACT 


Keywords: In-process detection of cutting forces, temperature, roughness, wear etc. during machining of 

Cutting force titanium alloys are very important. The Finite element (FE) analysis plays an important role in 

Cutting temperature monitoring and detection of machining responses. It offers a high accuracy in modeling of dry 

E aa a (FEM) cutting processes and its performance in modeling of cryogenic machining process is a matter of 

Titanium alloys machining interest. In this context, current investigation focuses on the dry turning and LN2/COz2 cooling 

In-process measurement assisted turning process of commonly used Ti6Al4V alloy. It is very useful material in the 
biomedical sector, and the simulation of cutting forces and cutting temperature via finite element 
method (FEM) has been performed. In addition, the simulation results are validated with 
experimental work. The results show that the deviations between FE modeling and experimental 
results for the cutting temperature are the average of 5.54%, 5.18% and 8.42% for the dry, LN2 
and CO% cooling conditions, respectively. On the other hand, the deviations from FE modeling 
and cutting force test results were 3.74%, 3.358%, and 3.03% under dry, LNz2 and COz cooling 
conditions, respectively. 
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1. Introduction 


The world of manufacturing especially machining sector is about to endure a seismic shift from conventional to advance processes. 
With this aim, Industry 4.0 supports the system that enables and reduces the cost of products [1]. Ultimately, obtaining a high quality 
product cannot be considered as a success alone [2]. However, achieving the most affordable production is among the main goal and 
expectation in the context of industrial revolution. To reduce the cost of the production process, the optimization of process parameters 
must be done well [3]. Optimization of machining parameters, cutting tool and cutting conditions reduce the cost [4]. Of course, in 
addition to all these, it is very common to demand a manufacturing carried out under conditions that may cause the least harm to the 
environment in the field of manufacturing [5-7]. 

These days application of sustainable cooling conditions i.e., dry, cryogenics especially with liquid nitrogen and carbon dioxide are 
used in machining of various materials [8]. These cooling conditions significantly improves the machining performance with less 
environmental effect [9,10]. Research studies have been carried out to solve the above expectations [11]. But every year, there are 
technological developments in cutting tools and cutting conditions used in the manufacturing process and machine tools [12,13]. In 
parallel with these developments, the necessity of continuing the research emerges. To reduce the machining time in the turning 
process [14], to obtain better quality surfaces [15], to reduce the cost and to minimize the environmental damage of the process 
[16,17], to improve and develop the cutting tools [18-20] used in turning have gained momentum in recent years. Nath and Kurfess 
[21] stated that cutting tools with different geometric features such as complex chip breaker form, tip angles of the cutting tip, cutting 
tool tip radius have been produced and widely used in the industry. However, by considering these geometric features in cutting tools 
as a whole with cutting parameters, it is necessary to focus on complex problems such as the effect of the chips formed during turning to 
reduce the damage to the cutting tool and the workpiece by breaking them in appropriate sizes, the workpiece surface quality and the 
effect on the cutting forces [22] and cutting temperatures [23] during cutting process. Thus, by determining suitable cutting condi- 
tions, it is aimed to contribute to obtaine products with minimal cost, time and optimal quality. With this context, there are lot of 
studies in the literature about finite element modelling of machining process by Papazoglou et al. [24], and especially dry conditions 
by Akgiin and Demir [25]. Zhou et al. [26] predicted the cutting forces in dry turning of Inconel 718 alloy with different tool nose 
radius and edge radius by finite element simulations. The authors verified the cutting forces by the FE model with the experimental 
studies in the range of 4% and 11%. Lotfi et al. [27] developed a friction model under 3D dry turning of AISI 304 stainless steel. The 
researchers verified that the created 3D friction model calculates the friction value depending on the cutting angle and force. Kumar 
et al. [28] created a 2D FE model to investigate the machinability performance of AISI 1045 steel via nanolayered coated WC-Co 
cutting insert in dry condition. It was found that the deviation between experimental and predicted force values varied from 1.8 to 
3.8% for Fx, which corresponds to the uncoated tool, but the values of coated tool is varied from 1.25 to 6%. Yang et al. [29] compared 
the surface roughness values obtained by experimentally and FE models in dry turning of Al6061 alloy by PCD cutting tool. They 
confirm that the good agreement between simulation and experimental responses has been observed, and the FE method can be used to 
simulate the surface roughness in related studies. Vijayaraghavan et al. [30] performed the FE simulations of dry turning of Inconel718 
alloys by data analysis method. The authors verified that the developed FE model is suitable to achieve a goal of sustainable 
manufacturing. Parida and Maity [31] carried out the hot turning of Inconel 625 alloy experimentally and supported it by FEM. There 
was a good relation and verification on chip formation and cutting forces between experimental values and FE results. Schindler et al. 
[32] predicted the thermal influences on Al2024 alloy by FE turning model. The researchers had a good accuracy between experi- 
mental and FE results based on the cutting forces and heat flows. Apart from above studies, a brief review about FEM studies for 
machining of difficult-to-cut materials are given in Table 1. 

As can be seen from studies in the literature, the FE studies were generally developed based on dry machining conditions. When 
cryogenic studies are performed, they are generally experimental, not finite element simulation type [41-45]. To the best of the 
authors’ knowledge, there are no or limited number of FE machining studies based on cryogenic cooling conditions in turning of 
Ti6Al4V alloys. Considering all these, the machining process is focused on the Ti6Al14V, by subjecting them to dry turning and LN» / 
CO» cooling turning process and the simulation of this process with the FEM has been performed. The findings and results are helpful 
for industrialists using such cooling technologies, and researchers conducting research on these issues. In addition, it has been 
observed that there are differences between the cooling technologies used in scientific studies and manufacturing industries. All studies 
examined in this section have shown that the finite element model simulations made by simulation softwares are possible to reach the 
results compatible with experimental studies. 


Table 1 

FEM studies for machining of difficult-to-cut materials. 
Material Machining Type Cutting Tool Machining Outputs References 
Ti6Al4V Turning —_— Cutting force, temperature Sadeghifar et al. [33] 
Inconel718 Turning Carbide Chip morphology, cutting force, temperature, contact length Razanica et al. [34] 
Ti6Al4V Turning Uncoated and TiAlN coated carbide Chip morphology, cutting force, stress, contact length Mishra et al. [35] 
Inconel718 Turning Tungsten carbide Cutting force, temperature, stress Parida et al. [36] 
Ti6Al4V Turning Uncoated coated carbide Chip morphology, cutting force, temperature Duet al. [37] 
Inconel718 Turning Tungsten carbide Chip morphology, cutting force, temperature, residual stress Silva et al. [38] 
Ti6Al4V Turning TiC coated PCBN Cutting force, temperature Veeranaath [39] 
Inconel718 Milling Uncoated ceramic SiAION Cutting force and Power Nieslony et al. [40] 
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2. Materials and methods 
2.1. Finite element method 


Finite element analysis is a numerical approximation method that allows complex structures to divided into particles and trans- 
formed into idealized structures and mathematically solved [46]. The finite element method includes systematic procedures of de- 
rivatives of finite elements that are sub-regions of approximation functions [47]. Numerical solutions for variables at special points are 
called nodes of each element and are applied to the calculation of finite solutions for the whole geometry [48]. Finite element analysis 
is useful in providing approaches for solving problems that are analytically difficult or impossible to solve [49]. Discretization of 
geometry is the first step in finite element analysis, and then the geometry is divided into finite elements. Next, different types of 
elements are used, such as triangular, quadrilateral, and shallow defined, with the element properties depending on what is required 
for the analysis. Stiffness matrix is used to define element properties [50]. This matrix expresses the force-displacement relationships 
of the element under load. Defining the stiffness matrix, the definition of the loading conditions such as pressure, force, and velocity 
throughout the boundary conditions at the specific node follows. Finally, the expression of the element, the applied loads and the 
boundary conditions are adapted to the formation of the matrix. The established equations are solved numerically with unknown 
values. Stress, strain and other properties of interest can be defined depending on the result of the displacement of the node. The 
equilibrium equations for general linear or nonlinear static problems of finite element analysis can be expressed as Eq. (1) [51]. 


{F} = [K]*{U} 1) 


Here, F is the vector matrix for the force at the element, K is the stiffness matrix and U is the vector of the determined node 
displacement. However, dynamical analysis requires a different form for dynamical problems independent of time as shown in Eq. (2) 
[52]. 


{FO} = Ku} + [m] 0) [me] fo} i 


Here “K” refers to the stiffness matrix, “U” to the node displacement vector refreshed with time change, “M 1” to the damping matrix, 
{uh to the node velocity (first derivative of the node displacement), ”M” to the mass matrix, and “Ü” to the node acceleration (the 


second derivative of the node displacement). In the formulation of finite elements; Lagrangian, Eularian, and Arbitrary Lagrangian- 
Eularian mesh structures are widely used. When these mesh structures are compared, it is seen that there are differences between 
them. In the Eulerian mesh structure, the finite element mesh remains constant during the material flow, only points are displaced. The 
benefit of the Eulerian mesh structure is that the shape of the elements does not change over time, so there is no distortion or distortion. 
However, it is imperative to consider the free surface of the chip as the starting shape and the chip formation process cannot be 
modeled. On the other hand, in the Lagrangian method, there is a movement in the roof. The net is fixed to the material and moves with 
the net material. However, it is a difficult situation to change the shape of the elements during the material flow. Therefore, the 
deformed web may need to be regenerated with the mesh. Thus, in the use of such mesh structures; the mesh needs to be re-meshed or 
adaptive meshing. Also, chip separation criteria are required for shaping chips. 

In the Arbitrary Lagrangian-Eularian (ALE) mesh structure, the points on the surface are not fixed neither to the material nor to the 
area, these points move arbitrarily. In other words, material flow is independent of points. According to Savidis et al. [53], ALE gives 
better results when compared with both mesh structures. ALE mesh algorithms are used in many softwares used for machining. For 
instance, it is used in softwares such as ABAQUS, DEFORM, Third Wave AdvantEdge, ANSYS LS-DYNA and gives better results. 
However, Third Wave AdvantEdge is an only machining simulation software although ABAQUS, DEFORM or LS-DYNA are the 
simulation softwares about general plastic deformation process. For that reason, Third Wave AdvantEdge was selected as software due 
to its being of specific FE simulation software that can easy to process. About this issue, Grzesik [54] indicated that the remeshing 
method used in some of these softwares increases the sensitivity and accuracy of the analysis. Some important boundary conditions are 
required to use the finite element method for analysis of chip formation in metal cutting. These boundary conditions help to define the 
limits and loading conditions in order to solve the problem. By applying these boundary conditions, determining the deformation rate 
of the material, and understanding how the workpiece material will behave during plastic deformation and it is possible to model the 
chip formation by selecting the appropriate modeling technique and modeling the material. To capture the dynamic behavior of the 
model, the strain rate and thermal effects of the constitutive equations should be selected. Stress in metals is changed as a function of 
strain, strain rate and temperature. Johnson-Cook constitutive equation is one of the first equations involving strain and thermal effects 
according to statement of Xu et al. [55]. Obviously, using these equations, material models such as rigid plastic, elastic-plastic, thermo- 
elastoplastic, thermo visco-plastic can be created [56]. The formulation of Johnson-Cook’s constitutive equation [57] used by Third 
Wave AdvantEdge is presented in Eq. (3). 


EP T-T, 
0 _ (A+B inl a 1- r \m 
o =(a+aey')(14+cn[ >] | (1-(F=)) 3) 


Here o° represents flow stress, æ plastic deformation, ¿ deformation rate, é refers to reference plastic deformation rate, T 
workpiece temperature, Tm workpiece melting temperature, Tr room temperature. The coefficient A indicates yield strength (MPa), B 
hardness modulus (MPa), C strain rate sensitivity coefficient, n hardness coefficient and m indicate thermal softening coefficient. 
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Moreover, for chip separation situation, a damage model is required. General JC damage model equation for equivalent fracture strain 
[58] is expressed in Eq. (4). 


r- (0,+vexos8)(1+040(£))(1+0,(E=2)) : 


In Eq. (4), 2 shows the stress triaxiality where P is the hydrostatic pressure that is the average of three-dimensional normal stresses. 


o 


z is the fracture strain during the chip separation which is a dynamic process of machining. The fracture constant of D,, D2, D3, D4, Ds 
in JC damage model equation shows the initial failure strain, exponential factor, triaxial factor, strain rate factor and temperature 
factor, respectively. The general geometric structure for finite element analysis with Third Wave AdvantEdge is presented in Fig. 1. As 
the turning method, 3D nose turning shown was preferred. 


2.2. Material details 


The numerical analysis of machining is directly related to the accurate adaptation of the material model into the software. 
Therefore, the success of modeling is possible by determining the material properties well and adapting them as input. During the finite 
element analysis, Ti6Al4V workpiece material having the dimensions of 1x1x3 mm size was selected. The JC parameters and fracture 
constants previously used by Xu et al. [55], chemical compositions by Gupta et al. [59] and the other mechanical features of Ti6A14V 
alloys used by Du et al. [37] in the study are given Tables 2-5, respectively. 

The cutting insert was selected as CVD (chemical vapour deposition) type TiC-Al203-TiN (each assumed as 1 um thickness [60]) 
coated carbide inserts with ISO coded of CNMG 120,408 as shown is Fig. 2. The properties of carbide cutting tools previously used by 
Du et al. [37] are provided in Table 6. 


2.3. Cooling systems 


Cryogenic is generally defined as the study of low temperature situations [61]. When analyzing a system operating under cryogenic 
conditions, it should be considered that the thermal features of materials change with temperature. Specific heat, thermal conductivity 
and especially heat transfer coefficient values for gases and liquids are important thermal properties for heat transfer calculations at 
cryogenic temperatures. The calculation of heat transfer per unit time (Q) is mentioned in Eq. (5) [62]. 

„dT 


Q = Cm (5) 


In Eq. (5), Cv, mand dT/dt specifies the specific heat, the mass and the rate of temperature decrease per unit time, respectively. By 


OET 


substituting “m = A*6*p” in Eq. (6); 
dT 
— CR AK SK ne 
Q = C,*A*5*p at (6) 
The heat flux of density “q” equals to “Q/A”. So, Eq. (6) take its the form of Eq. (7) by Newton’s heat transfer laws; 


= *C age (7) 
q=p iy dt 


From Eq. (7), dT/dt is known as the ratio of the temperature differences like the cooling rate (T2-T1)/(t2-t1). 


q 


h= 
T, Ty 


(8) 


In Eq. (8), “h” is the heat transfer coefficient (W/m?K), Tw is the surface temperature and Tris the cryogen temperature of LN2 and 
CO as —196.5 °C and —78.5 °C, respectively. Moreover, the heat transfer coefficients for LN and CO, are adapted as 309 W/mK and 


Cutting speed (V) 


Fig. 1. General geometric structure for finite element (3D nose turning). 
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Table 2 
JC parameters of Ti6Al14V alloys [55] 
A(MPa) B(MPa) n C m 
782.7 498.4 0.28 0.028 1 
Table 3 
JC fracture contents of Ti6Al4V alloys [55] 
D1 D2 D3 D4 D5 
—0.09 0.25 —0.5 0.014 3.87 
Table 4 
The chemical content of Ti6Al4V alloys [59] 
Ti Al vV Fe c oO 
89.9% 5.8% 4.0% <0.25% <0.3% <0.2% 
Table 5 
The other mechanical and thermal properties of Ti6A14V alloys [37] 
E(GPa) Tm(°C) a(10/°C) k(W/mK) v p(kg/m?) ¢p(J/kg°C) 
108 1650 9.1 6.8 0.33 4430 526.3 


Select Insert Parameter 
1.Insert Shape 2.Clearance Angle 3. Tolerance 4.InsertType  5.InsertSize 6. Insert Thickness 7.Nose Radius 
N v A A 00 00 08 v 
Inset Shape 
Side Rake Angle -6 (deg) [_]Wiper Geometry 
Wiper Radius 3 (mm) 
Back Rake Angle 0 (deg) z 
80° Wiper Offset 0.2 (mm) 
Lead Angle 0 (deg) `na” Wiper Offset 
Edge Radius 0.04 (mm) e Y 
Tool Width 1 (mm) 
“x _ Wiper Radius 
ATs 
Advanced Options Cancel Help { S No% Radius 
Fig. 2. Tool insert parameters of CNMG 120408. 
Table 6 
Mechanical and thermal properties of the carbide tool (Copyrights reserved) [37] 
E(GPa) k(W/mK) a(10°°/°C) v p(kg/m?) ¢p(J/kg°C) 
640 50 4.5 0.22 11,900 220 


165 W/m?K [63], respectively as in Fig. 3. 


2.4. Meshing 


The mesh structure created in the finite element analysis directly affects the analysis results. The main goal is to select a mesh 
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Coolant Options 


Coolant 
© Constant Properties 
Heat Transfer Coefficient [h] {Wé(m2.K) } 309 
Coolant Temperature {degC} -196.5 Coolant Area Immersed 


© Temperature Dependent Properties 


Input Table 


Coolant Area Option 


®© 


© Exclude Tip Vicinity 

Distance from Cutting Edge Tip (R) {mm} 0 ô S 3 
I 
I 


Heat Flux (q) 


© Focused Location {mm} 
| Workpiece 
x J0 Y |0 Z 10 Lal (0,0) E 
Jet Radius (R) {mm} 0 In the Immersed Option, the heat flux BC is applied on non-contacting 


surfaces. 


Fig. 3. Cooling options for LN2 and CO3 cryogens. 


structure that can represent the geometry as much as possible. The density of FEs is an essential measure used to confirm the accuracy 
of the analysis (element type and morphology also affect the accuracy of the analysis). Assuming the model has no singularity region, 
the higher density mesh produces outcomes with high precision. Though, the element mesh of the model is very dense and requires a 
lot of computer memory and long execution time. This drawback often occurs under multiple repeating conditions, especially for 
nonlinear and transitional analysis. Comparing the findings to test data or theoretical values is one approach to assess the quality of a 
finite element mesh. Unfortunately, in the early stages of research, test data and theoretical conclusions are frequently unavailable. 
Therefore, more different tools are required to assess mesh quality. The most important one is the meshing sensitivity analysis. For the 
current study, the minimum element size for the workpiece and the cutting tool was optimized as 0.15 mm by the mesh sensitivity 
analysis with re-meshing method shown in Table 7. The simulations were performed with a computer having the RAM of 16 GB, the 
client cache of 6 GB and the processor speed of 2.4 GHz. All the 8 cores of the computer were used for parallel processing to solve FE 
analysis faster. It means that 8 FE simulations can be performed at the same time in same solving period. Firstly, the element sizes of 1, 
0.5 and 0.25 mm were used in FE analysis and then the results were obtained as 320.1, 298.4, 285.3 N with a computing time of 0.7, 1.5 
and 3.3 h, respectively. This showed that the element size should be more decreased due to higher values of change in cutting forces in 
small computing time differences. Afterwards, the element size was decreased to 0.1 mm step by step. The change in cutting force is so 
small (ignorable) but the computing time increased about 70%. For that reason, the element size was determined as 0.15 mm (Fig. 4) 
with a computing time of 6.1 h for a simple test simulation. 


2.5. Friction model 


In FE simulation of the turning process, Özel [64] suggested that the Coulomb friction model may be quite efficient on the tool flank 
surface. The average friction coefficient between tool and chip in orthogonal turning was computed from by measuring cutting force as 
presented in Eq. (9) from a simple turning test for each condition. As a result of the calculation, the average friction coefficient (u) was 
found to be 0.46, 0.42, 0.37 for the dry, LN2 and CO% conditions, respectively. In Eq. (9), Fe and Ft are the main cutting and tangential 
forces obtained experimentally with a dynamometer while a is the rake angle. For instance, average friction coefficient (u) is calculated 
as 0.46 when Fc, Ft and « are 270.5 N, 159 N and —6°, respectively. Tangential forces were measured but not recorded in previous 
study of Gupta et al. [63] 


Table 7 
The mesh sensitivity analysis for a simple analysis. 
Element size (mm) Cutting force (N) Simulation solving time (h) 
1 320.1 0.7 
0.5 298.4 1.5 
0.25 285.3 3.3 
0.2 276.6 4.8 
0.15 270.5 6.1 
0.1 269.5 10.2 
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Minimum 
element size 


Maximum 
element size 


Fig. 4. The minimum element size for a simple test simulation. 


_ F, + F.*tana 
— F, — F*tana 


(9) 


3. Experimental procedure for verification 


In this work, the same experimental conditions of Gupta et al. [63] were adopted for performing the simulation. The machining of 
commercially used Ti6A14V alloy was performed with the CVD type TiCn-Al203-TiN coated carbide with ISO coded of CNMG 120,408 
rhombic shaped inserts. The work piece dimensions were the diameter of 30 mm and the length of 100 mm. The initial layer of 0.5 mm 
was removed in order to maintain uniformity of workpiece. The side cutting edge angle also known as principle cutting angle used was 
90° having 0.8 mm nose radius. Further, the experiments were performed under same conditions with FE simulations. For this purpose, 
the CNC lathe machine with control system of Siemens were used. The specifications of CNC center and parameters used for cryogenic 
cooling system by Gupta et al. [63] are given in Table 8. 

The levels of the feed rate and the cutting speed under dry, LN2 and CO% were selected based on the previous work of Gupta et al. 
[63]. As mentioned earlier, the same machining conditions were adopted and the experimental design is shown in Annexure 1. The 
cutting depth of 0.5 mm and the machining time of 60 s were used throughout the experiments. These parameters were purely selected 
by following tool manufacturer guide and preliminary studies. The details of machining parameters were listed in Table 9. Then, the 
cutting temperature and cutting forces were evaluated with the help of Fluke make thermal camera and Kistler made 9257A dyna- 
mometer. All these measurements were in online mode and the dynoware package was utilized to extract the cutting forces from the 
system. The complete process methodology is given in Fig. 5. 


4. Results and discussions 
4.1. FE simulation results 
The simulation study performed with Thirdwave Advantage software using the process parameters and boundary conditions 


detailed are shown in above sections, 3D turning was analyzed and a new finite element model was developed. Fig. 6 shows a sample 
test of finite element simulations for the feed rate of 0.1 mm/rev and the cutting speed of 50 m/min. 


Table 8 
Technical parameters of CNC turning center and cooling conditions [63] (Copyrights reserved). 
Machine Specifications 
CNC lathe machine Minimum and maximum spindle speed 200-3000 RPM with 6KW motor. 
Cooling condition Flow rate of coolant: 0.35 L / min, pressure: 4 bar, dimensions of nozzle: 1 mm brass nozzle diameter, 2 nozzles with position at rake 
details face 
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Table 9 
The machining factors and their levels [63] (Copyrights reserved). 
Factors Level 
1 2 3 
Cutting speed, V (m/min) 50-75—100 
Feed rate, f (mm/rev) 0.1-0.15-0.2 
Cooling conditions Dry - LN2 - CO2 


4.2. Main cutting forces 


In metal cutting processes, cutting force is the force necessary to remove material from a resistive body. The principle cutting force 
and its parts are affected by a variety of variables. Cutting parameters, workpiece material, cutting tool properties and cutting 
environment are the most basic factors known to be effective on cutting force according to Yang etal. [65]. Knowing the cutting force is 
helpful when it comes to sizing the cutting tool, estimating the power required for cutting, optimizing the cutting parameters, etc. 
Considering all this, as a result, minimal cutting force is often desirable in machining operations. In current section, the influence of 
cutting parameters and cutting environment on the main cutting force has been determined via finite element analysis and the results 
are presented graphically in Fig. 7. In parallel with the increasing cutting speed, an increase in cutting force was observed in all three 
feed rates and cutting environments. For example, when the calculation is made according to each increasing level of cutting speed, the 
increase in cutting force was between 16.56% and 23.70% at 0.1 mm/rev feed, between 18.24% and 28.17% at 0.15 mm/rev feed, and 
between 14.23% and 26.36% at 0.20 feed mm/rev. In many previous studies, unlike this study, with increased cutting speed, Pathak 
et al. [66] and Şeker et al. [67] proved that there is a reduction in cutting force in machining of Al-(1-2) Fe-1 V-1Si alloy and St44 steel, 
respectively. Ravi and Pradeep Kumar [68] stated that the reason for this situation is the thermal softening effect of the work material 
because of heat generated by rise in cutting speed in machining of titanium alloys. However, in current work, the force ascends with 
elevated cutting speed. Generally, when the workpiece materials are subjected to deformation, strain hardening occurs in the material 
in addition to the thermal softening event. An increase in cutting force can be observed if the strain hardening effect exceeds the 
thermal softening effect. Therefore, in this study, the rise in force versus cutting speed was associated with strain hardening as a similar 
observation was also done by Hou et al. [69] and Shi et al. [70]. As can be seen from the Fig. 7, a rise in feed triggered an increase in 
cutting force. Depending on the rise in feed, the increase in cutting force ranged from 6.72% to 23.11%. However, when compared to 
the influence of cutting speed, the feed had lower impact on cutting force. This phenomenon may be associated with higher feed rate, 
which increases the contact area at the tool-chip interface as more material is removed per unit time mentioned by Kim et al. [71]. 
When Fig. 7 is examined, it can be determined that the cutting environment offering the lowest cutting force in each cutting parameter 
is LN2 cooling followed by CO2 and dry environment. Compared to dry cutting, 21.01% to 34.95% less cutting force can be achieved 
with LN3 assisted cutting. In addition, thanks to CO3 assisted cutting, this improvement can reach up to 14.41%. The improved strength 
and hardness of the tool material, decreased tool wear, and eliminated adhesion between the tool-chip and tool-workpiece interfaces 
gained by reducing the cutting temperature may be result in lower cutting forces achieved by LN3 and CO% assisted machining. Since 
LN; has a lower boiling temperature than COg, it prevents the cutting tool from overheating during machining and offers more effective 
cooling ability as thought by Pusavec et al. [72]. Also, Hong et al. [73] indicated that the outstanding performance of LN2 cooling may 
be that the liquid/gas buffer across the contact surface absorbs the heat generated and provides a lubricating effect that evaporates 
quickly and reduces friction. 


4.3. Cutting temperature 


Cutting temperatures are highly dependent on machining parameters, workpiece material and tool material characteristics. In 
particular, the characteristics of the material being machined are of vital importance in determining the range of cutting temperatures. 
For example, Singh et al. [74] demonstrated that titanium alloys, which have much lower thermal conductivity than steels, lead to 
higher cutting temperatures during machining. Because Ti6Al4V has a limited thermal conductivity, roughly 80% of the heat 
generated stays in the cutting zone. On the other hand, during the machining of some materials like aluminum, Pramanik [75] declared 
that more heat is thrown out along with the chips. In the current section, the effect of cutting parameters and cutting environment on 
the main cutting force has been determined via finite element analysis and the outcomes are provided graphically in Fig. 8. The upward 
effect of cutting speed on temperature can be noticed from Fig. 8. The degree of this effect varied between 23.61% and 81.25%, 
depending on the feed rate and cutting environment. In general, a higher cutting speed accelerates the build-up of cutting temperature. 
It was resulted of more heat produced per unit time and less time available for heat evacuation from the region [76]. Feed rate was 
found to have a lesser effect on the cutting temperature. That is, the increase in temperature was up to 16.89% in response to an 
increase in feed from 0.1 to 0.2 mm/rev. This is most likely due to the increased feed causing more friction between the removed 
material and the cutting tool, which increases the energy in the system and hence raises the cutting temperature. This finding is 
consistent with prior research of He et al. [77]. With machining under LNg, the cutting temperature is drastically reduced compared to 
machining under dry and CO2 machining. This phenomenon is due to the fact that when LN3 is injected to the tool-chip contact region, 
it absorbs heat and evaporates rapidly from area, causing the cutting temperature to drop. LN2 lowered the temperature from 58.43% 
to 73.53% compared to dry, and from 40.32% to 57.14% compared to CO2. Another observation is the slowing down of the LN2 
temperature reduction rate with increasing cutting speed and feed. This is evidence that the cooling efficiency of LN2 decreases at 
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elevated cutting speeds and feed rates. In addition, since LN3 has a lower boiling temperature than CO2, Pusavec et al. [72] stated that 
it prevents the cutting tool from overheating during machining and offers more effective cooling ability. 

Moreover, to verify the data in Fig. 8, the steps of a simple finite element simulation for the evaluation of the maximum cutting 
temperature are given in Figs. 9, 10 and 11 for the dry condition, LN and CO3 cooling conditions, respectively. 

As shown in Figs. 9-11, there is a maximum temperature on the contact surface of the cutting tool. For this reason, the heat 
generated between the tool and the chip affects the performance of the tool and the surface quality of the machined workpiece. The 
friction that occurs in the slip zone also affects the mechanical properties of the machined material and causes tool wear, and higher 
cutting forces due to the heat generated in the tool, chip-tool and workpiece contact area. Moreover, the relation between cutting 
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Fig. 6. Sample finite element simulations for the feed rate of 0.1 mm/rev and the cutting speed of 50 m/min. 


forces and cutting temperature well explained in Fig. 12. 


4.4. Verification of FE analysis results by turning experiments 


In Fig. 12, the comparison of the cutting temperature and cutting forces obtained as a consequence of the experimental study and FE 
analysis shown graphically. The results obtained are very close to each other and the curves obtained from the graph quite consistent 
with each other. This phenomenon shows that the experimental results verified the FE model results in high accuracy. 

Based on Fig. 12a-c-e, the deviations between FEM and experimental results for the cutting temperature are the average of 5.54%, 
5.18% and 8.42% for the dry, LNz and CO% cooling conditions, respectively. On the other hand, the deviations from FEM and cutting 
force test results in Fig. 12b-d-f averaged 3.74%, 3.35%, and 3.03% under dry, LN2 and COz2 cooling conditions, respectively. In the 
literature, most of the FE machining studies shows that the deviation less than about 10% is thought as the acceptable range between 
the finite element model results and experimental validation as mentioned before by Özel et al. [78]. After the verification of the FE 
model, the possible wear region of the cutting inserts was evaluated for each cutting conditions. Fig. 13 shows the possible wear region 
of the inserts evaluated by FE model and the actual wear region of the cutting inserts after experimental turning previously used by 
Gupta et al. [63]. The experimental and FE images of the cutting inserts are also so similar with respect to tool wear, and this phe- 
nomena also verified the FE model which is developed for the cryogenic turning process of Ti6Al4V alloy. 
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Fig. 7. Finite element results for the main cutting force at the feed rate of a) 0.1 mm/rev, b) 0.15 mm/rev, c) 0.2 mm/rev. 


5. Comparison of the results 


The most common reason about the deviations between the FE model and experimental results are the meshing accuracy and 
friction model since these parameters are just assumed in the FE model. The other possible reason for this is considered to be that the 
Johnson Cook parameters adapted for the workpiece material model are obtained from the literature. However, the mechanical 
properties of the workpiece material used in the present study may have different from the same materials in the literature due to 
different microstructures caused by different manufacturing methods or heat treatments as thought by Guo et al. [79]. Mechanical 
properties are mainly due to interatomic bond strengths. However, the internal structure (microstructure) of the material also has an 
effect on it. In this way, it becomes possible to obtain different mechanical properties, namely tensile, compressive strength, hardness, 
toughness in the same material by changing the microstructure. Moreover, there may be different tensile strength values of the same 
material under the same condition of the tensile tests. For that reason, it is generally taken the average of the 3 ~ 5 tensile strength 
values from the identical situation of tests [80]. On the other hand, one of the other methods differing the microstructure of the same 
material is precipitation hardening. Abdel-Salam et al. [81] defined the precipitation hardening as to ensure that the second phase, 
which is present in less amount in the material and it is precipitated in the form of particles in the main phase, thereby increasing the 
strength of the material. According to Fu et al. [82], this method is one of the most important methods used to increase the strength of 
materials and is generally used in non-ferrous metal alloys (Al, Ti, Mg). The main reason for the increase in strength in the precipitation 
hardening process is explained by the limitation of the dislocation movements of the precipitates formed as a result of the precipitation 
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Fig. 8. Finite element results for the cutting temperature at the feed rate of a) 0.1 mm/rev, b) 0.15 mm/rev, c) 0.2 mm/rev. 
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Fig. 10. Finite element steps for the cutting temperature analysis of the insert at the feed rate of 0.1 mm/rev and the cutting speed of 50 m/min in 
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Fig. 11. Finite element steps for the cutting temperature analysis of the insert at the feed rate of 0.1 mm/rev and the cutting speed of 50 m/min in 
CO» conditions. 


of the supersaturated melt. However, the adapted material in the finite element software may behave just as with the certain adapted 
parameters, namely JC parameters, damage constants, tensile parameters, elastic modulus, etc. The special conditions into the 
structure of the material cannot be defined for the material in the software. For that reason, the average deviations of 6 ~ 7% may be 
attributed to these circumstances. Korkmaz et al. [57] investigated the cutting force in dry turning of Nimonic 80 alloy as both 
experimentally and FE model and had about 6% deviations between the results. The researchers commented on the reason for the 
deviation that it may be attributed to both the microstructure and hardness of the workpiece material. They also stated that the de- 
viations may be assigned to the experimental ploughing effect due to smaller depth of cut than tool nose radius. However, it is not 
possible to model this effect in FE software. For the cryogenic conditions, the heat transfer coefficients are assumed as constant in the 
FE software. Wang et al. [83] thought that the cutting and cooling process is a thermodynamically process, and thus it should be 
actually modelled dynamically or temperature dependent type but then Zimmerschied and Isermann [84] demonstrated that it will 
take so much time to solve a simple FE simulation in cryogenic conditions. That is why in the current FE study, the heat transfer 
coefficients and cooling temperatures of LN2 and CO2 cryogens were adapted as constant. This issue may be attributed to the deviations 
between FE model and experimental results in cryogenic conditions. 


6. Conclusions 


This study focuses on the dry turning and LN2/COz cooling assisted turning process of Ti6A14V and the simulation of this process by 
modeling it with the finite element method. The aim of this study is to to verify the FEM simulation study with the experimental turning 
process and to show the estimation of the results that can be obtained from the machining process. The prominent results of this study 


are summarized as follows: 


e The force increased with increasing cutting speed. The increase in force versus cutting speed was associated with strain hardening. 
Depending on the increasing cutting speed, the increase in cutting force was between 16.56% and 23.70% at 0.1 mm/rev feed, 
between 18.24% and 28.17% at 0.15 mm/rev feed, and between 14.23% and 26.36% at 0.20 feed mm/rev. On the other hand, an 
increase in feed triggered an increase in cutting force. Depending on the increase in feed, the increase in cutting force ranged from 


6.72% to 23.11%. 


13 


M. Kumar Gupta et al. Mechanical Systems and Signal Processing 169 (2022) 108772 


== Deviations —Ħ— Simulation —Ħ— Experimental = Deviations —Ħ— Simulation —®— Experimental 

450 

420 
© z 
g 390" 
z E a : 
£ pi £ s60 £ 
ž & & be 
g E È 330 f3 
a bn © 
£ & 
= G] 

E 300 
o 2 

270 

o 3 6 9 12 15 18 21 24 27 0 3 6 9 12 15 18 21 24 27 


Experimental runs Experimental runs 


(a) (b) 


[C] Deviations —Ħ— Simulation —"— Experimental [Deviations —*— Simulation —®— Experimental 
330 10 360 


14 


12 300 8 330 
= 
A 7 
10 
270 6 300 
8 5 
m 240 a 270 
m, 
s n ‘ 3 
210 240 
2 
=" “N | = 
1 — 1 
o 3 6 9 12 0 3 6 9 12 15 18 21 %4 27 


% 
Cutting temperature, °C 


Deviations 
Deviations 


a 


Main cutting forces, N 


4 
2 
1 
o o 
15 18 21 24 27 
Experimental runs Experimental runs 
(c) (d) 
Deviations —Ħ— Simulation —Ħ— Experimental |__| Deviations —Ħ— Simulation —Ħ— Experimental 
is qS S SIS 7 390 
164 64 
+ 180 + 360 
144 t 
Lo 5 z 
124 £ P L 330 $ 
a | ie = 2 ¢ 
5 150 = E a4 £ 
£104 P] Ě be 
f= g E = 
b p = . = 
Z g z -300 .E 
á °] E 8% z 
j} 120 & G 
64 E- „m E 
= 24 H270 3 
44 5 z 
24 gs a | 240 
o — T T o 
o 3 6 9 12 15 18 21 2% 27 v 3 6 9 12 15 18 2 2% 27 
Experimental runs Experimental runs 


(e) (9 


Fig. 12. Cutting temperature and cutting forces comparison of the FEM and experimental results; a-b) Dry, c-d) LN2, e-f) CO2. 
e It can be determined that the cutting environment offering the lowest cutting force in each cutting parameter is LN2 cooling fol- 


lowed by CO% and dry environment. Compared to dry cutting, 21.01% to 34.95% less cutting force can be achieved with LN2 
assisted cutting. In addition, thanks to COz assisted cutting, this improvement can reach up to 14.41%. 
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Fig. 13. The possible wear region of the inserts evaluated by FE model and the actual wear region of the cutting inserts after experimental turning 
for; Dry, b)LN2, and c)CO2 (Adopted and Modified: Copyrights reserved) [63]. 


e The upward effect of cutting speed on temperature was observed. The degree of this effect varied between 23.61% and 81.25%, 
depending on the feed rate and cutting environment. Feed rate was found to have a lesser effect on the cutting temperature. That is, 
the increase in temperature was up to 16.89% in response to an increase in feed from 0.1 to 0.2 mm/rev. 

e With machining under LNg, the cutting temperature is drastically reduced compared to machining under dry and CO2 machining. 
LN2 lowered the temperature from 58.43% to 73.53% compared to dry, and from 40.32% to 57.14% compared to CO2. Another 
observation is the slowing down of the LN2 temperature reduction rate with increasing cutting speed and feed. This is evidence that 
the cooling efficiency of LN2 decreases at high cutting speeds and feed rates. 
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e The experimental results verified the FE model results with high accuracy. The deviations between FE modeling results and 
experimental results for the cutting temperature are the average of 5.54%, 5.18% and 8.42% for the dry, LN2 and CO% cooling 
conditions, respectively. On the other hand, the deviations from FE modeling results and cutting force test results were 3.74%, 
3.35%, and 3.03% under dry, LN2 and CO2 cooling conditions, respectively. 

e Although the deviation between the experimental results and the FE model results is within acceptable limits, it is not possible to 
mention about a complete agreement. Since the heat transfer coefficients and cooling temperatures of LN2 and CO2 cryogens were 
adapted as constant, this issue may be attributed to the deviations between FE model results and experimental results in cryogenic 
conditions 

e This research will help to research & development centers of the machining industries , especially working on improvement of 
cooling technologies in the machining of biomedical materials. 
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